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Abstract 

Thermal model fit indicates early chemical freeze-out of multi-strange hadrons with small collec- 
tive velocities at 200AGeV Au-Au collisions at RHIC. In this work, we present our recent results 
by SPheRIO hydrodynamical calculations inspired by this picture. In our model, multi-strange 
hadrons go through chemical freeze-out when the system reaches some temperature close to the 
phase transition, stopping to make inelastic collisions, and their abundances are therefore deter- 
mined only by partonic EOS. At a lower temperature thermal freeze-out takes place where elastic 
collisions are brought to a halt. We calculate the spectra for various hadrons at different centrality 
windows, with chemical and thermal freeze-out temperature being fit as a function of centrality. As 
it is shown, the result provides a reasonable panoramic description of the spectra of identified par- 
ticles. Chemical freeze-out gives good correction of the multiplicity of certain species of particles, 
especially for multi-strange hadrons. 

PACS numbers: 
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I. INTRODUCTION 



In heavy-ion collision a hot and dense matter is formed which eventually evolves into a 
state of freely streaming particles. Strong collective flow patterns measured at RHIC sug- 
gest that the hydrodynamical model is well justified during the intermediate stages of the 
reaction. The model provides a useful tool for drawing information of space-time evolution 
of the matter, from the time when local equilibrium is attained until the process of hadron 
decoupling, called freeze-out. The simplest scenario to treat freeze-out is to adopt chem- 
ical/thermal freeze-out temperature |l|, |2], it is widely used to give qualitative estimation 
of overall properties of system. However, it is worth noting that conclusions drawn about 
characteristic of chemical/thermal freeze-out temperature are essentialy based on statis- 
tics and hydrodynamics-motivated model fittings, (i) chemical freeze-out: Statistics and 
hydrodynamics-motivated models give us the characteristic temperature of T c h in relativis- 



tic heavy-ion collisions. It is determinec 



the experimental data 



from fits of particle yield ratios for hadrons to 



10 



11 



12, [l3, Q, 15, l3, 17| . The temperature 



indicates when the hadron abundances become fixed. Fits are done for various centrality 

3, 14, 15, 16, 17]. For 
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windows, and reveal that T c h is almost independent of centrality 
200AGeV Au-Au collisions, the value of chemical freeze-out temperature is around 160MeV, 
close to the phase boundary. However, it is usually assumed strangeness is not fully equi- 



librated, thus a strangeness saturation factor 75 is introduced. Un 



ike chemical freeze-out 



temperature, 75 varies with collision system and centrality 15j, ll7|]. The compilation of 
chemical freeze-out parameters (T ch , fif h ) at various collision energies seem to be aligned on 
one curve in the T — fi B plane. And can be parametrized by the average energy per particle 
< E > / < N > ~1 GeV|6j. (ii) thermal (kinetic) freeze-out: The thermal (kinetic) freeze- 
out temperature T t h is obtained from the slope of transverse momentum distribution with 
some radial flow profile. On the thermal freeze-out surface, elastic rescattering processes 
cease and hadrons start to escape freely. Statistical model fit shows that temperature at 
the kinetic freeze-out depends on the centrality Ij], 0, [3, Q, 20, 21]. It is generally under- 
stood that the system undergoes first the chemical freeze-out where the observed particle 
ratios are fixed and next the thermal freezeout where the shape of the transverse-momentum 
distribution is fixed. Thus T c h>T t h is expected. 

It is worthwhile to check to what extent chemical/thermal freeze-out picture as well as 
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above conclusions can be confirmed by the more realistic and accurate full hydrodynamic 
model, where the space-time evolution of fluid is taken into account rigorously. So in this 
work we explore the scenario of chemical freeze-out by calculating transverse momentum 
spectra of 200A GeV Au-Au collisions. A full three-dimensional hydrodynamical model with 
chemical freeze-out mechanism is employed. The results are compared with experimental 
data for all centrality windows. 

The paper is organized as follows. In the next section, we briefly review the hydrody- 
namic model we used and explain how chemical freeze-out is implemented. Our results and 
discussions are presented in the last section. 



II. SPHERIO WITH CHEMICAL FREEZE-OUT 



The hydrodynamical model we adopted is based on smoothed particle hydrody- 
namic(SPH) algorithm 23l. |24j, In this model, the matter flow is parametrized in terms 
of discrete Lagrangian coordinates, of the so-called SPH particles. As a result, the hydro- 
dynamic equations are reduced to a system of coupled ordinary differential equations. The 
main advantage of the model is that it can tackle the problem with highly asymmetrical 
configurations, as is concerned in relativistic high-energy nucleus-nucleus collisions. 

The code which implements the entropy representation of the SPH model for relativistic 
high energy collisions, and which has been developed within the Sao Paulo - Rio de Janeiro 
Collaboration, is called SPheRIO. It has been successfully used to investigate the effects 



of the initial-condition fluctuations and adopting t 



re continuous emission scenario for the 



description of decoupling process 
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27|. 



Usually, highly symmetric and smooth initial conditions is employed in hydrodynamic 
approach. However, owing to the fact that in heavy ion collisions the system is small, 
fluctuations are not negligible. To take this into account, we use an event simulator, NEXUS, 
to generate the initial condition of SPheRIO. This provides us possibility to study the 
collision on fluctuating event-by-event basis. 

In the previous applications of SPheRIO, the strangeness had been neglected. Here, 
SPheRIO has been further improved in order to consider strangeness conservation and to 
adopt the scenario of chemical freeze-out. The strangeness conservation is implemented by 
explicitly incorporating strangeness chemical potential into the code, and correspondingly a 
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new set of equation of state(EOS) has been built and utilized. We use the hadronic resonance 
model with finite volume correction to describe the matter on the hadronic side, where the 



main part of observed resonances in Particle Data Tables [221] has been included. For quark 
gluon plasma (QGP) phase, the ideal gas model is adopted. As a good approximation, we 
assume local strangeness neutrality throughout the hydro evolution. 

As it was shown |28j]. hydrodynamic model without chemical freeze-out gives good de- 
scription of transverse momentum spectra of light hadrons. While it reproduces the shape 
of spectra for hyperons and anti-proton, there are visible discrepancies in the multiplicities. 
To compensate this, we implement chemical freeze-out for strange hadrons such as A, S, 
Q and and anti-protons. At chemical freeze-out temperature T c h, these particles cease 
to have inelastic collisions and therefore their abundances are fixed. We assume complete 
chemical equilibrium on the surface of chemical freeze-out, thus no strangeness saturation 
factor 75 is introduced in our treatment. With the system being cooled and rarefied further, 
thermal freeze-out occurs at a lower temperature T t h, where the distribution function and 
abundances of the rest of hadrons are determined. In our calculation, T c h and T t h serve as 
adjustable parameters as functions of centrality according to experimental data. 

III. RESULTS AND DISCUSSIONS 

The thermal freeze-out temperature is tuned as a function of centrality, to reproduce the 
shape of transverse momentum distribution of all charged particles. Furthermore, as in the 



previous works 



23 



24 



25 



26(, a rescaling factor is introduced to fix the pseudo-rapidity 
distribution for all charged particles. 

In Fig.l, we computed the pseudo-rapidity distribution for all charged particles for 
Au+Au collisions at 200A GeV. It is depicted for six different centrality windows. The 



experimental data are from PHOBOS Collaboration 29]. 

The experimental transverse-momentum-distribution data for all charged particles is pre- 
sented in Fig. 2, the experimental data are reproduced with a choice of freeze-out temperature 
from T t h = 135 MeV for most central collisions to T t ^ = 150 MeV for most peripheral ones. 
We use linear interpolation to determine the thermal freeze-out temperatures of intermedi- 



ate windows. The experimental data are from STAR Collaboration 
200A GeV, with 77 = 0. 



301 ]. taken in Au+Au at 
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With the parameters chosen as explainded above, we calculate in the following the spectra 
for various hadrons. In Fig. 3-8 we show the transverse- momentum spectra of pions, protons 



and kaons for different centra 



ity windows at mid rapidity, as well as experimental data 



from BRAHMS Collaboration 31]. The spectra of A, H and Q are depicted in Fig. 9-13, 
together with data from STAR Collaboration 32]. We use dotted lines to represent the 
results obtained without incorporating chemical freeze-out, and solid lines for those with 
chemical freeze-out switched on. The interpolated values of thermal and chemical freeze-out 
temperatures are noted in the figures. 

It is observed that the present hydrodynamic model gives good description of the exper- 
imental transverse-momentum spectra for pions, kaons and protons even without turning 
the chemical freeze-out on for these particles. Without the chemical freeze-out incorporated, 
although it gives the correct slopes for the spectra of anti-protons and strange hyperons A, 
S and Q, the disagreement comes from the multiplicities of the spectra. By introducing 
chemical freeze-out, those results are improved significantly. While the spectra for pions, 
kaons and protons almost remain the same, (this is not shown in the figures.) it provides a 
good fit of anti-proton and strange hyperons A, 5 and Q. The chemical freeze-out increases 
the multiplicities of strange hadrons with respect to the ones in the original model owing 
to higher chemical freeze-out temperature, as compared with the thermal freeze out one. 
Meanwhile, the slopes of the spectra remain unchanged. It was indicated experimentally in 
Ref. 33j], as one goes to large rapidity region where the baryon density differs sizably from 
zero, the outcome of chemical freeze-out would be much more significant. Further work on 
this topic is under progress. 

It is worth noting that the effects of chemical freeze-out have been discussed by several 



authors. In ref 



34 



35 



36 



371 ]. the early CFO is studied in terms of hybrid model where the 



hydodynamic evolution is complemented with a hadronic cascade model where the value of 
chemical freeze-out temperature is not studied explicitly. Other hydrodynamic calculations 
have been carried out either only in the transverse plane assuming Bjorken's scaling 38( or 
with the baryon chemical potential taken to be zero [391] . and only some of the centrality 
windows are covered. 
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FIG. 1: Pseudo-rapidity distributions of all charged particles for Au+Au collisions at 200A GeV. 
The data are from PHOBOS Collaboration [29f]. 
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□ STAR Data, ii=0 all charged particle 
SPheRIO, T fit: (MeV) 




FIG. 2: Transverse-momentum distributions of all charged particles for Au+Au collisions at 200A 
GeV in the pseudo-rapidity interval —1.0 < r/ < 1.0. The data are from STAR Collaboration 30]. 
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FIG. 3: Transverse-momentum distributions of tt for Au+Au collision at 200 A GeV in the pseudo- 
rapidity interval —1.0 < rj < 1.0. The data are from BRAHMS Collaboration [3 1 ] . 
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FIG. 4: Transverse-momentum distributions of n + for Au+Au collisions at 200A GeV in the 
pseudo-rapidity interval —1.0 < r\ < 1.0. The data are from BRAHMS Collaboration [3 1 ] . 
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FIG. 5: Transverse- momentum distributions of K for Au+Au collisions at 200A GeV in the 
pseudo-rapidity interval —1.0 < rj < 1.0. The data are from BRAHMS Collaboration [3 1 ] . 
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FIG. 6: Transverse- momentum distributions of K + for Au+Au collisions at 200A GeV in the 
pseudo-rapidity interval —1.0 < r\ < 1.0. The data are from BRAHMS Collaboration [3 1 ] . 
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FIG. 7: Transverse- momentum distributions of protons for Au+Au collisions at 200 A GeV in the 
pseudo-rapidity interval —1.0 < rj < 1.0. The data are from BRAHMS Collaboration [3 1 ] . 
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FIG. 8: Transverse-momentum distributions of anti-protons for Au+Au collisions at 200 A GeV in 
the pseudo-rapidity interval —1.0 < r\ < 1.0. The data are from BRAHMS Collaboration [3 1 ] . 
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FIG. 9: Transverse-momentum distributions of A for Au+Au collisions at 200A GeV in the pseudo- 
rapidity interval —1.0 < rj < 1.0. The data are from STAR Collaboration 32| . The dotted lines 



indicate the results obtained without chemical freeze-out, and solid lines for those with chemical 
freeze-out incorporated. 
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FIG. 10: Transverse-momentum distributions of A for Au+Au collisions at 200A GeV in the pseudo- 
rapidity interval —1.0 < rj < 1.0. The data are from STAR Collaboration 32| . The dotted lines 



indicate the results obtained without chemical freeze-out, and solid lines for those with chemical 
freeze-out incorporated. 




P T (GeV) 

FIG. 11: Transverse- momentum distributions of H + for Au+Au collisions at 200 A GeV in the 
pseudo-rapidity interval —1.0 < r] < 1.0. The data are from STAR Collaboration 32j . The dot- 



ted lines indicate the results obtained without chemical freeze-out, and solid lines for those with 
chemical freeze-out incorporated. 



12 



> 

CD 



T3 



10" 
10' 
10 a 
10" 3 
10" 4 
10"' 
10' 
10"' 
10' 





O s- iilAK Data, n=0 




T h =135(MeV), chemical freeze-out off 


] 


T ch /T th fit: (MeV) 




1 45/1 35(0-5%) , 1 48/1 38{ 1 0-20%), 1 51/141 (20-40%), 




~*-^eL^ 1 55.6/1 45.6(40-60%),1 60/1 50(60-80%) 








- - - ~-x S ^V~_ 






: o 

' - v ° 




s — 1 

\ 




T 

1 


D 

1 1 1 1 1 1 ■ 1 ■ 1 ■ 



1.5 2.0 2.5 

P T (GeV) 



FIG. 12: Transverse-momentum distributions of H for Au+Au collisions at 200A GeV in the 
pseudo-rapidity interval —1.0 < 77 < 1.0. The data are from STAR Collaboration 32j . The dot- 



ted lines indicate the results obtained without chemical freeze-out, and solid lines for those with 
chemical freeze-out incorporated. 
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FIG. 13: Transverse-momentum distributions of f2 + f2 + for Au+Au collisions at 200A GeV in 
the pseudo-rapidity interval —1.0 < r/ < 1.0. The data are from STAR Collaboration 32]. The 



dotted lines indicate the results obtained without chemical freeze-out, and solid lines for those with 
chemical freeze-out incorporated. 
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